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All galaxies without a radio-loud AGN follow a tight correlation between their global FIR and 
. radio synchrotron luminosities, which is believed to be ultimately the result of the formation of 

massive stars. Two face-on colliding pairs of galaxies, UGC 129 14/5 and UGC 813/6 deviate from 
this correlation and show an excess of radio emission which in both cases originates to a large 
extent in a gas bridge connecting the two galactic disks. The radio synchrotron emission expected 
from the bridge region is calculated, assuming that the kinetic energy liberated in the predom- 

*n ' 

■ inantly gas dynamic interaction of the respective interstellar media (ISM) has produced shock 

waves that efficiently accelerate nuclei and electrons to relativistic energies. A simple model for 
the acceleration of relativistic particles in these shocks is presented together with a calculation of 
the resulting radio emission, its spectral index and the expected high-energy gamma-ray emission. 
This process is not related to star formation. It is found that the nonthermal energy produced in 
the collision is large enough to explain the radio emission from the bridge between the two galax- 
ies. The calculated spectral index at the present time also agrees with the observed value. The 
expected gamma-ray emission, on the other hand, is too low by a factor of several to be detectable 
even with foreseeable instruments like CTA. 
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1. Introduction 

The universal, tight correlation between the spatially integrated far-infrared (FIR) luminosities 
and the monochromatic radio continuum emissions from late type galaxies without a bright, radio- 
loud AGN (de Jong et al. 1985, Helou et al. 1985) has been known for a long time. It is believed to 
be ultimately a result of the formation of massive stars and also holds for most interacting galaxies. 
It is therefore quite unusual that Condon et al. (1993), hereafter referred to as CHS, and Condon et 
al. (2002) found two clear exceptions to this rule. This concerns two face-on colliding spiral galaxy 
systems where at the present time, presumably some 30 Myr after the interaction, the respective 
pairs of galaxy disks are again well separated from each other optically, but are connected by a 
radio continuum-bright "bridge" of gas, suggested to be stripped from the interpenetrating disks. 
The two systems show overall a significant excess, by a factor of about two, of the radio continuum 
emission relative to the FIR-radio continuum ratio expected from the FIR-radio correlation for 
single galaxies. CHS interpreted the finding as the result of general electron escape from galaxies, 
in this case into the bridge connecting the pair. 

In the present paper the dynamical effects of such galaxy-galaxy collisions on the interstellar 
gas are investigated. It is argued that the interstellar media of the respective galaxies will undergo 
a largely gas dynamic interaction, where the low-density parts exchange momentum and energy 
through the formation of large-scale shock waves in the supersonic collision. A simple model 
for the acceleration of relativistic particles is presented and the synchrotron emission from the 
relativistic electron component is calculated, as well as the expected gamma-ray emission from 
relativistic nuclei and electrons. The model is shown to be able to explain the radio continuum 
emission observed from the bridge between the galaxies. For details, see Lisenfeld & Volk (2010), 
hereafter referred to as LV. 

2. Further properties of the colliding systems 

The spectral index of the radio emission between 1.49 and 4.9 GHz steepens gradually from 
the stellar disks with values of 0.7-0.8 to values of 1.3-1.4 in the middle of the connecting gas 
bridge. This steep spectral index will be argued to be indicative of dominant synchrotron and 
inverse Compton losses suffered by the relativistic electrons. Apart from the sychrotron emitting 
electrons the bridge in the system UGC 12914/5 contains also large amounts of atomic (CHS) 
and molecular (Braine et al. 2003) gas, where practically no star formation is taking place. This 
is interpreted as an essentially complete hydrodynamic removal of the more diffuse atomic and 
molecular gas from the galaxies. Most likely, only the dense cloud cores - capable of forming stars 
- have remained within the stellar disks together with the stars. The other system, UGC 813/6, was 
described in the later paper by Condon et al. (2002). Since it is very similar, the discussion here 
will be limited to UGC 12914/5. 

3. Acceleration model 

In a face-on collision the stellar disks interpenetrate each other without being too much altered. 
However, the diffuse gas and part of the gas clouds interact hydrodynamically and exchange energy 
and momentum. 
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If one assumes that half of the gas, which is now present in the bridge, was previously in one 
galaxy, and the other half in the other galaxy, then the total energy liberated in a fully inelastic 
interaction is the kinetic energy of the gas mass: 



where M gas is the total gas mass in the bridge and v co u denotes the velocity difference of the gas 
at collision. The factor 1/2 converts this velocity to the velocity difference in the center of mass 
system (assuming that both gas disks are equally massive). The relative velocity between the 
galaxies, v co u ~ 600 km/s, in the case of UGC 12914/5 has been derived by CHS from a analysis 
of the HI line and of the galaxy masses. M gas and E^ n are estimated as (4 — 11) x 10 9 M and 
therefore (0.49 - 1.4) x 10 58 erg, respectively. 

At the collision of the two ISM, a tangential discontinuity will form and two strong shocks 
will propagate in opposite directions with velocities v s h oc k, communicating the interaction to larger 
and larger fractions of the colliding interstellar gas masses. The space in between these shocks is 
filled with post-shock gas. Fig. 1 shows the idealized picture of this interaction, the basis of the 
present model, in the reference frame of the motion normal to this tangential discontinuity, which 
is also the center of mass system. In this frame, the post-shock normal velocity v post vanishes and 
the preshock normal velocity of the gas is v pre = 5V co u. The contact discontinuity is stationary and 
situated in the middle between the galaxies (at x = in Fig. 1). In a strong, but approximately 
unmodified, adiabatic gas shock, the normal component of the velocity difference between the 
shock in the pre- and postshock gas follows the relation: 4 x (v post — v s h oc k) = (vpre — Vshock)- With 
v pos t = 0, this yields v s h oc k = gVcoii- The shock velocity V s , relative to the unperturbed ISM gas, 
is then V s = v s h 0C k + 2 v coii = § v co ii = 400 km/sec. From Fig. 1 in CHS one can estimate that the 
shocks are at present close to the galaxy disks so that practically the entire bridge is expected to be 
filled with post-shock gas. 

The Mach number of the shocks produced in this collision is like in middle-aged supernova 
remnants in the Sedov phase. Roughly speaking, the particle acceleration efficiency of such shocks 
will therefore be similar to that of a supernova remnant, i.e. of the order of 10 — 30% (e.g. Berezhko 
& Volk 1997). This is a basic assumption for the present paper. 

The source function Q(E), i.e. the number of relativistic particles produced by the shock per 
energy and time interval is given by: Q(E) = 2f acc (E) v s h oc k A, where A is the area covered by 
the shock (roughly the area of the galaxy disks). The factor 2 is due to the fact that the shocks 
propagate into two, opposite directions. f acc (E) is the downstream, uniform number of relativistic 
particles of rest mass m, produced per volume and energy interval at the shocks: 



Here, y is the spectral index of the differential relativistic particle source spectrum, taken to be y = 
2.1 (Drury et al. 1994, Berezhko & Volk 1997). For the nuclear particles, essentially protons, the 
constant fo can be determined by requiring that the total energy converted into relativistic particles 
during the entire duration of the interaction, T, is equal to E acc = x E^ B , where = 0.1 — 0.3: 




(3.1) 




(3.2) 
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Figure 1: Schematic illustration of the shock in the bridge region, in the center of mass system. 
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for 7 > 2. Here, E m [ n is the minimum energy of relativistic protons of mass m p accelerated, taken 
as E mi n = 1 GeV (Drury et al. 1994). This gives 



/o = ^, acC ' (3.3) 



gacc (7-2) 

2rAv shoc k (m p c 2 ) 

In the Galactic cosmic rays, at a given energy, the number of relativistic electrons is about 1% 
of that of the protons at GeV energies, and the source spectra for electrons and protons are probably 
similar (Miiller 2001). Assuming that this electron-to-proton ratio is also representative for the 
source spectra in the galaxies considered here, implying that the electron and proton source spectra 
are the same, for the source function of the relativistic electrons Q e (E) = Q{E) x 0.01 is adopted, 
which means that their source distribution is / acc e = /o,e (E/(m p c 2 )) ~ 7 , with /o,e = fo x 0.01. 

4. Relativistic electron density and synchrotron emission 

In order to calculate the synchrotron emission from the bridge the time-dependent propagation 
equation for the electron particle density f e (t,x,E) is solved. Due to the spatial symmetry of the 
situation, a one-dimensional approximation is appropriate, where x is the coordinate in the direction 
along which the galaxies separate (see Fig. 1). In addition, the diffusion of relativistic electrons is 
neglected because the typical spatial scales which are relevant on the time-scales discussed here, 
3 x 10 7 yr, are < 1 kpc, whereas the width of the bridge is 10 kpc. Then one can write: 

df,(t£,E) = q ^ E) + d_ f bE i MttXtE) \ (41) 
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where 

q e (t,x,E) = -j-(8(x-v shock t) + 8(x + v shock t)) (4.2) 

is the local source strength (in units of relativistic electrons produced per energy interval per time 
and per volume). This source strength describes two shocks that start at t = at x = and propagate 
into opposite directions with velocity v s h 0C k- Eq.(4.1) takes into account the electron acceleration in 
the shocks together with the radiative energy losses of these CR electrons due to Inverse Compton 
and synchrotron losses; up to the present epoch, adiabatic losses can be disregarded. Then: 



where Oj is approximated by the Thompson scattering cross section, B the magnetic field strength, 
£/b its energy density and U m d is the energy density of the radiation field. Neglecting adiabatic 
losses is consistent with the fact that the distribution of the synchrotron spectral index a(l .49, 4.86), 
between 1.49 and 4.86 GHz, has an integrated value a(1.49,4.86) > 1 in the bridge region (CHS), 
indicative of dominant Inverse Compton and synchrotron losses. 

The magnetic field can be estimated from the minimum energy requirement to be about 7 jlG 
in the bridge of UGC 12914/5 (CHS), which implies U B = 1.2 eV cm" 3 . The total radiation 
field density from (in essentially equal amounts) the blue magnitude, the FIR, and the Cosmic 
Microwave Background (CMB) amounts to U m d ~ 0.76 eV cm~ 3 (LV). 
The solution of eq. (4.1) is: 

A( £ 'W-') = A.(^)" , {'-«('-w)r- <4 " 4) 

for t > |x|/v shock and t - |x|/v shoc k < Ti oss 

and 

f e (E, \x\,t) = 0, for t < |x|/v sh0 ck and t - |x|/v shoc k > Ti oss , (4.5) 

with Ti oss = (Eb)~ l being the life-time of a relativistic electron against radiative energy losses. 
This expression can also be integrated over the volume of the bridge to obtain the total number of 
relativistic electrons in the bridge F s (E,t). The result is: 

F e (E,t) = 2A/ 0e ( — ^-j] Vshock^loss— ^ " {1 - (1 - ^ Aioss) 7-1 } , for ? < Ti oss (4.6) 
\ m p c z / 7 — 1 

and 

F e (E,t) = 2A/ ,e j Vshock^loss r, for t > Z loss (4.7) 

The synchrotron spectrum is obtained by convolving f e (E,x,t), respectively F e (E,t), with the syn- 
chrotron emission spectrum of a single electron. 

This solution allows a theoretical prediction for the total radio flux density specifically at 
1.49 GHz and of the radio spectral index between 1.49 and 4.86 GHz. Both values are in satisfac- 
tory agreement with the observations (LV). Clearly also the radio spectral index increases towards 
the center of the bridge. The model calculations therefore show that acceleration by large-scale 
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shocks caused by the galaxy-galaxy interaction is indeed able to explain the radio emission from 
the bridge both in intensity and spectral morphology. 

The two interacting galaxy pairs studied here, are possibly examples of what might have hap- 
pened much more frequently at early stages of structure formation, when primordial galaxies had 
already developed magnetic fields as a consequence of early star formation, but when they were 
still likely to interact strongly with neighboring structures of a similar character. 

5. High-energy gamma-ray emission 

Although not the primary topic of this paper, it is clear that the interaction of galaxies con- 
sidered here will also lead to the acceleration of gamma-ray producing very high-energy particles, 
both nuclei and electrons, in the form of the distribution / acc , cf. eq. (3.2). The visibility of the 
acceleration process also in high-energy gamma rays would be an independent argument for the 
model presented. In the following a rough estimate will be given. 

The shock system which characterizes the face-on interaction of the two spiral galaxies can 
to first approximation be considered as plane parallel, with constant speed shocks; the energetic 
particles remain confined in the growing interaction region until the interaction is completed. As- 
suming the diffusion coefficient to be as low as the Bohm diffusion coefficient (e.g Kang 2007) 
one can approximately calculate the maximum proton energy achieved at the present epoch as 
£max ~ 3.7 x 10 16 eV. This very high energy is the result of the long life-time (T = 2.8 x 10 7 yrs) 
of the shocks, inspite of their comparatively moderate speed of pa 400 km sec -1 . The situation 
is rather different for the energetic electrons because of their radiative losses. Even disregarding 
any magnetic field amplification at the shocks due to the accelerating particles, the maximum elec- 
tron energy is only ~ 10 TeV. For the IC gamma-ray emission in the sub-TeV region mainly the 
CMB counts. In contrast to the great majority of the radio electrons, however, all gamma-ray emit- 
ting electrons are "old" as a result of postshock radiative losses. As a consequence their energy 
spectrum is softened with a correspondingly reduced IC gamma-ray emission. 

When the shocks have gone through the interacting ISM of the two galaxies, which is the case 
at about the present epoch, the kinetic energy E^ n = (0.49 to 1.4) x 10 58 erg has been transformed 
into thermal and nonthermal particle energy. This corresponds to E acc ~ 10 57 erg, predominantly in 
relativistic nuclei, which is roughly 10 7 times more energy than available from a single supernova 
remnant. 

The bridge volume V is reasonably estimated as V TtR 3 , where the disk radius R 10 kpc 
(CHS). This implies a mean gas density of 0.05 to 0.14 cm~ 3 , which is an order of magnitude 
smaller than the typical density in the plane of spiral galaxies and more than three orders of mag- 
nitude smaller than the density in the starburst nucleus of e.g. the galaxy NGC 253. 

An analytical estimate for the integral hadronic gamma-ray emission, from ^"-production by 
collisions of energetic protons with gas nuclei and subsequent decay into two gamma rays, is given 
in Eq. (9) of Drury et al. (1994) for gamma energies E large compared to 100 MeV: 

F ( > E,, = T) »9 X ,0-0 (^) (JL) - 1 X p h „«„„ S cm - s - 
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Inserting the values E^ n = 10 58 erg, d = 61 Mpc, = 0.1, and n = 0.05 to 0.14 cm 3 results in 

(_1) £ 2^X>^) ^ p tQ 6) x 10 -15 erg cm -2 (5 1} 

approximately independent of energy up to about 10 15 eV. 

For the IC emission of the electrons at high gamma-ray energies the distribution F e (E, t) for t > 
Ti oss in eq. 4.7 is relevant. For y « 2 the spectral energy flux density is roughly energy independent 
with a value ~ 7 x 10~ 15 erg cm~ 2 s _1 . 

Thus the IC gamma-ray energy flux is of the same order as the hadronic gamma-ray energy 
flux in the region of energy overlap, even for the radiatively cooled electrons. This is a consequence 
of the low mean gas density in the bridge. On the other hand, the lowest TeV photon flux from an 
astrophysical source detected until now was F[> 220 GeV) = 5.5 x 10 -13 photons cm~ 2 s _1 . For 
a flat spectral energy density this corresponds to 8.8 x 10~ 13 erg cm~ 2 s . The measurement was 
made with the H.E.S.S. telescope system for the nearby starburst galaxy NGC253 (Acero 2009). 
Taking this result as a yard stick, the expected hadronic flux from UGC 12914/5 at gamma-ray 
energies above 1 TeV is still two orders of magnitude below this minimum flux. The expected 
minimum detectable energy flux for the future Cherenkov Telescope Array (CTA) is as low as 
«5x 10~ 14 erg cm~ 2 s _1 in the TeV region, and at least one order of magnitude higher at 50 
GeV (CTA 2010). Therefore the gamma-ray flux is also below the detection capabilities of CTA. 
To this extent the more optimistic expectation by LV is corrected here. This flux is also expected 
to be too low for the detection capabilities with the LAT instrument on Fermi at lower gamma- 
ray energies. The same source at the distance of NGC 253 would be detectable even for present 
Northern Hemisphere ground-based TeV gamma-ray instruments like VERITAS and MAGIC, and 
at GeV energies for Fermi. Ultimately the reason for the low gamma-ray fluxes from these distant 
interacting systems is their low gas density and their comparatively high age, despite the large 
overall nonthermal energy they contain. 
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